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We investigate theoretically shot noise in an artificial Single Molecule Magnet based upon a CdTe 
quantum dot doped with a single S=5/2 Mn spin and gated in the hole sector. Mn-hole exchange 
anisotropy is shown to lead to profound consequences on noise, like super-Poissonian behaviour. We 
report on a novel effect, similar to the Dicke effect in Quantum Optics, that allows to separately 
extract the hole and Mn spin relaxation times using frequency-resolved shot noise measurements. 
We expect that our findings may have further relevance to experiments in other S > 1/2 systems 
including transport through Mn 12 molecules and STM spectroscopy of magnetic atoms. 
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Single Molecule Magnets (SMMs) combine properties 
of a magnet with those of a nanostructure. This makes 
SMMs a useful platform to merge concepts and applica- 
tions of spintronics and nanoelcctronics 1]. In a tran- 
sistor setup, SMMs behave as Quantum Dots (QDs). 
This allows detailed Coulomb Blockade (CB) level spec- 
troscopy of prototypical examples, like Mn^fl, 0| or 
endofullerene N@Ceo The experimental extraction 
of their intrinsic properties is, however, very challenging 
because little is known about the precise role of various 
sample-preparation steps like the attachment of leads. 
An alternative route is to study artificial counterparts to 
a SMM, like a CdTe QD doped with a single Mn+ 2 ion 
with spin S — 5/2 (Fig. 1). This proposal is based on 
recent experiments where these QDs (without contacts) 
are optically probed @, @, 0, H, Q ■ When doped with a 
single hole, the system behaves as a SMM with magneti- 
zation steps and hysteresis (lo| . Furthermore, its trans- 
port properties depend on the quantum state of the Mn 
spin,giving rise to remarkable phenomena like hysteretic 
CB |H|. 

While dc transport entails considerable information, a 
complete understanding requires to go beyond and study 
shot noise (HI, which we address here. Our results show 
that exchange anisotropy leads to super-Poissonionan 
shot noise due to lack of relaxation of the Mn spin. Inter- 
estingly, a novel effect at finite frequencies, similar to the 
Dicke effect in Quantum Optics [131 ] , allows to separately 
measure the hole and Mn spin relaxation times (Fig. 4). 

Model. — The minimal model of an artificial SMM 

reads H® f D f = J2a £ h d U<? + H Tff ■ Thc first term de- 
scribes confined holes within the QD, the second term 
describing the effective hole-Mn anisotropic exchange is 
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FIG. 1: (Color online) Schematics of the artificial SMM based 
on a CdTe QD doped with a single Mn +2 ion and coupled to 
metallic reservoirs (tunnelling rates Fl and Tjj and chemical 
potentials Hl/r — Vo is the gate voltage. 



6|3/2, —1/2) and | U) = a|3/2, -3/2) + 6|3/2, +1/2). 
H®R is obtained from the full problem H = J2 n £ n,d\ l d n + 

JhSh(rM)M after projecting onto the | fh), I JJ-ft) sub- 
space [llf . d^ n creates a confined hole in the spin-orbital 
ip n (f), described by a six band Konh-Luttinger Hamilto- 
nian, and Shivn) is the hole spin density evaluated at 
the Mn location. J = Jh\^{fu)?, Jh ^ 60eVA 3 being 
the hole-Mn exchange coupling constant of CdTe. Due to 
strong spin-orbit coupling, exchange is highly anisotropic 
and spin flips are suppressed (JJ./j \S^\ fT/i) ~ unless 
there is some heavy hole-light hole mixing due to asym- 
metry b ^ 0, a = J x /J ^ [HE!. 

Quantum Master Equation (QME). — The full 
Hamiltonian including the coupling to reservoirs reads 



H = H 
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The Hamiltonian 



of each hole reservoir is H" es = J2k a a £k c 

whereas H% = J2k a a^a c k + H.c. describes tun- 
neling. The QME for the reduced density matrix, pit), 
is obtained after applying a Born-Markov approxima- 
tion 11 1 with respect to H^ R , pit) = Cp(t). Thc 
dissipative dynamics is governed by a Liouvillcan su- 
peroperator C which contains the forward/backwards 



(±) transition rates Nzfll 



N 



spin operators and r^ :± are Pauli matrices operating in -EW^i) J2 a | (iV | | iV =F 1) | 2 between eigenstates of H^R 



thc hole lowest energy doublet | ft h ) = o|3/2, +3/2) 
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The couplings T^n are assumed to be constant and 
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N=0 

- 5/2, 0),| + 5/2,0) 

- 3/2, 0),| +3/2,0) 
I 1 to n\ l -X. i /9 n\ 


E = 
£ = 


N=l, antiferromagnetic (AF) 
-5/2,1^,1 + 5/2,^) 
|-3/2,fr fe >,| + 3/2,^) 

1/9 -fV, \ 4-1/9 lit \ 


E = e h - 5/4 J 
E = e h - 3/4 J 
W, — fi —1/4 7" 


N=l, ferromagnetic (FM) 

-l/2,J| h ),| + l/2,^) 
|- 3/2, U), \+3/2,th) 
\-5/2,U),\ + 5/2,tih) 


£ = s A + 1/4,7 
B = e h + 3/4J 
E = £ll + 5/4J 



TABLE I: States for the pure Ising case (Eq. JTJ with J ± = 0) . 
N=0: the six M z projections of the Mn spin are degenerate. 
N=l: the presence of a single hole breaks the degeneracy. 



/+ {fa = 1 — fa) 1S the Fermi function In the 

strong CB regime we only consider states with N=0,1, 
holes (Table I). The steady state p stat , is obtained as 
p(t) = Cp stat = 0, such that C has a zero eigenvalue with 
right eigenvector |0}} = p stat . The corresponding left 
eigenvector is |0)) = i such that <<0|0>> = Tr[lp stat ] = 1 
[II, and averages read (A) = Tr{Ap stat } = ((0|i|0)}. 

Shot Noise. — The shot noise spectrum S(uj) — 
2 f™ x d T e lu}T {{AI{t+T), A7(t)}) is defined as the Fourier 
transform of the irreducible (cumulant) correlation func- 
tion ((I(t + r)J(t))) = (AI(t + r)AI(t)>, where AI(t) = 
I(t) — (I) measures deviations away from the steady 
state current (I). With the help of the MacDon- 
alds formula [la ]. S(u>) can be rewritten as S(u>) — 
u> J °° dt sin (uit)-^((n 2 (t))). {(n 2 )) is the variance of the 
number n of holes being transferred in, say, the right 
lead in the time interval t. This process is stochastic and 
is governed by the probability distribution P n (t) which 
gives all moments (or cumulants) like, in particular, 
((n 2 )). In this context, P n (t) is expressed in terms of the 
so-called n-resolved density operator, p^ n \t), as P n {t) = 



Tr sys {p n (t)} 0. 



' sys 
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(t) fulfills a generalized QME 
•+Jn-1) _,_ /--„(n+l) 



pv> = (C - £+ - Cr)p^ + Cy^-V + C R p [n+1 K The 
superoperator £^ contains forward/backwards rates. Us- 
ing the above method, the noise at the right barrier reads: 

Srr(w) = ((0\J\0)) - ((0\1R(lu)1\0)) - (u -to), (2) 

X = — is the current superoperator (such that 
(I) = ((0|T|0))) whereas J = + jC^ is a superoper- 
ator describing self-correlations at the barrier [l9j. The 
pseudoinverse operator is defined as R(oj) = Q i( jL c Q, 
where Q = 1 — |0))(0| projects out the null space de- 
fined by |0)} = p stat [181 ] . Taking into account the 
left barrier Sll(m>) and the noise coming from charge 
accumulation [2(1 . described by the cross-correlations 
Slr(u),Srl{w) [21(, the total noise (Tj, = Tr) reads 
S(lo) = \{Sll{uj) + Srr(uj) + Slr(u) + Srl(u)}- 

Results for J± — 0. — The shot noise properties in 
the J± = case are investigated in Fig. 2. The top 



panels show results for the ui — shot noise (Fig. 2a) 
and Fano factor F = 5(0)/2(7) (Fig. 2b) as a func- 
tion of V and Vq. At finite bias voltages, 5(0) presents 
steps which depend on Vq. Between plateaus, 5(0) 
changes at values of V where the differential conduc- 
tance Gdc{V) = d(I)/dV has maxima (not shown) such 
that fluctuations are enhanced. In the limit V — > 0, 
our calculation recovers the fluctuation-dissipation the- 
orem 5(0) = 4kTGdc(0), which relates thermal fluctu- 
ations (Jonhson-Nyquist noise) with the linear conduc- 
tance. In this linear response regime, both shot noise 
and linear conductance exhibit a three-peak structure 
in regions of gate voltage corresponding to charge de- 
generacy between the N = and N = 1 sectors (from 
Vq = —3 to Vq = —5 in the figure). This is in stark 
contrast with a normal QD which would exhibit instead 
one single CB peak. Key for an understanding of this 
unusual CB is the fact that M z does not relax during 
transport 11]: In the absence of holes, the spin of the 
Mn ion is free and therefore all the six projections M z 
are degenerate (Table I, top). Sweeping the gate voltage 
towards the charge degeneracy region, the only allowed 
transitions are those which conserve M z . This results in 
three possible charge degeneracy points corresponding to 
the transition A r = 0<^>A r = lat different values of 
M z . For example, the first CB peak (Vq = —3.75) cor- 
responds to charge degeneracy between | — 5/2,0) and 
| - 5/2,117,) (or |5/2,0) and |5/2,^ h ». The charge de- 
generacy condition for states with either M z — ±3/2 or 
M z = ±1/2 occurs at higher \V G \ (Table I, middle) and, 
therefore CB is spin-dependent. By further decreasing 
Vq, we obtain a very small, but finite, three-peak struc- 
ture which can be attributed to thermal excitations to 
excited states in the ferromagnetic (FM) sector (Table 

1, bottom). The above physical picture leads to strong 
spin-dependent deviations from Poissonian noise (Fano 
Factor, Fig. 2b). For \Vq\ < 5, the noise remains sub- 
Poissonian due to spin-dependent CB. For |Vg| > 5, the 
competing dynamics between slow (filled states in the AF 
sector) and fast channels (FM sector) results in bunching 
and the noise becomes super- Poissonian [22] ■ As V in- 
creases, FM excited states enter within the bias window 
and then the noise becomes sub-Poissonian. In this case, 
bunching of holes at low voltages is replaced by direct 
hole spin- flips due to tunneling to/from the QD. The re- 
sults for finite to are shown in the bottom panels of Fig. 

2. At Vq — —3.75 (Fig. 2c), shot noise is frequency- 
independent up to V ~ 3. In this voltage region, only 
states corresponding to the AF sector come into play In- 
terestingly, at V > 3 a resonance around u> = develops. 
The appearance of this resonance can be understood as 
an enhancement of fluctuations due to tunneling through 
both AF and FM channels which are now within the bias 
window. For example, starting from the state | — 1/2, fa), 
one hole can tunnel out of the QD leaving the Mn spin 
in the state M z = —1/2. A second hole can now tunnel 
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FIG. 2: (Color online) Noise of an artificial SMM for a = 0. Top panels: Zero-frequency Shot noise (a) Fano factor (b). 
Bottom panels: Finite frequency Fano factor, 5(w)/2(J) for Vg ~ —3.75(c) and Vg ~ —5.5(d). Parameters: eh = 5, J = 1, 
T L = F R = 0.01 and T = 0.05 (All energies in meV) 
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FIG. 3: (Color online) Zero- frequency Fano factors for 
a) Vg = —3.75, b) Vg = —5.5 (the rest of parameters are 
the same as in Fig. 2). Bottom panels: blow up of low Fano 
factor regions of the corresponding top panels. 



onto the state | — 1/2, JJ-^), which is energetically avail- 
able, resulting in an effective spin relaxation for the hole 
(M z does not change). This mechanism gives rise to a 
resonance in S{oj) whose width is given by the inverse 
hole spin- flip time 1/T[ L k, T (other intrinsic hole spin- 
flip mechanisms, not included here, would also contribute 



to this width). Further reduction of Vg allows both kind 
of channels to participate in transport even at low volt- 
ages and the resonance is present for all V (Fig. 2d), in 
good agreement with our previous interpretation. 

Role of spin- flip terms. J± ^ 0. — If valence band 
mixing is nonzero, |±3/2) heavy holes couple with |=pl/2) 
light holes ( J± ^ in Eq. |T]) and b ^ in the expressions 
for | f|-/j) and | JJ-/i)). This mixing allows simultaneous spin 
flips between the hole and the Mn spins [1, [23| which re- 
sults in split states which are bonding and antibonding 
combinations of \M Z — —1/2,^) and \M Z = +1/2,^). 
This splitting can be extracted directly from the d(I)/dV 
curves (not shown here). Fig. 3 illustrates the dramatic 
changes that mixing induces on noise. As a = J±/J 
increases, the Fano factor can reach values F >> 1 in 
voltage regions (both gate and bias) where the noise is 
sub-Poissonian for a = (Fig. 3a). Interestingly, the 
largest Fano factors are obtained for small values of a: 
Most of the time, transport events conserve M z . How- 
ever, spin-flips mediated by hole transport are small but 
nonzero. As a result, periods of small current, followed 
by larger currents, are possible leading to huge Fano fac- 
tors. We emphasize that super-Poissonian behavior is 
related here with Mn spin-flips. This is in contrast with 
the a = case where only hole spins can flip. Although 
we do not seek here a detailed explanation of the vari- 
ous F >> 1 regions for different values of Vg and V, we 
mention in passing that each particular feature in Figs. 
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3 can be associated with different transitions between 
tt^O states. 

S(lo) and the Dicke effect. — After many transport 
cycles occur, the Mn spin relaxes completely in some typ- 
ical time scale Tf^ » T/ 1 : Starting from an empty QD 
with the Mn in the spin state M z , a hole tunnels and ex- 
changes one unit of spin with the Mn. After some time 
t m 1/r, the hole tunnels out of the QD such that the 
Mn spin is now M 2 ±l. Note that it takes many tun- 
neling events to completely relax the Mn spin, such that 
Ti >> Ti- This separation of time scales is reflected 
in S(lo), which, remarkably, consists of a narrow peak 
on top of broad resonance (Fig. 4a). The explanation 
of this feature is a process analogous to the Dicke effect 
in Quantum Optics ID, 2^ ■ the coupling of both relax- 
ation channels (fast spin hole relaxation and slow Mn 
spin relaxation) leads to a splitting into two combined 
decay channels for the whole system. The width of the 
superradiant channel (broad resonance) allows to extract 
the hole spin relaxation time as l/T[ l . More importantly, 
the subradiant channel (narrow resonance) can be used 
to extract the intrinsic Mn spin-flip time as 1 /Tf 1 in an 
independent manner. When a — 0, only hole spin relax- 
ation is possible (Fig. 4a, solid line). As a changes, a 
extremely narrow resonance develops on top of the broad 
one. Furthermore, the difference in magnitude between 
the noise values at uj = (see previous section) and 
1/Ti > lu > l/Ti ' 1 may facilitate the experimental ver- 
ification of this effect. The broad resonances for a = 
and a ^ are very similar, as one expects from our inter- 
pretation in terms of hole relaxation. When one of the 
spin relaxation channels is not available, the Dicke ef- 
fect should disappear. We have made systematic checks, 
as a function of V, Vq and a, which corroborate this 
and demonstrate the robustness of the effect. Fig. 2c, 
gives a clear demonstration of the absence of Dicke effect 
when a = 0. The opposite case, namely when hole spin 
relaxation starts to be inefficient, should result in the 
disappearance of the broad feature in the noise (Fig. 4b). 

Concluding remarks and experimental consequences. — 
In summary, our calculations show that Mn-hole ex- 
change anisotropy has profound consequences on shot 
noise, like super-Poissonionan behavior due to incom- 
plete relaxation of the Mn spin. The main result of this 
Letter is a novel Dicke effect in the shot noise spectrum 
that can be used to separately measure the spin relax- 
ation times of, both, holes and Mn. Finally, most of the 
physics captured by our model is inherent to large spin 
S > 1/2 systems with anisotropy. We therefore expect 
that shot noise measurements in other SMMs, like in the 
experiments of Refs. 0, 0, Q , m ay reveal the effects de- 
scribed here. STM spectroscopy of magnetic atoms such 
as S = 5/2 Mn on C1L2N [25| is one further experimental 
example where our findings may be relevant. 

We thank Joaquin Fernandez-Rossier for his input 
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FIG. 4: (Color online) a) Dicke effect in the shot noise spec- 
trum S(ui) of an artificial SMM. The width of the superradi- 
ant channel (broad resonance) allows to extract the spin hole 
relaxation time as l/7\ h whereas the subradiant channel (nar- 
row resonance) gives the Mn spin relaxation time as l/T^ • 
b) At lower voltages, hole spin relaxation mediated by charge 
fluctuations is inefficient resulting in a single narrow peak due 
to Mn spin relaxation. For clarity, we do not show the largest 
values at u) = 0, see main text. 
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